Electron Bunching in Stacks of Coupled Quantum Dots 
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We study the transport properties of two double quantum dots in a parallel arrangement at 
temperatures of a few Kelvin. Thereby, we show that decoherence entailed by the substrate phonons 
affects the shot noise. For asymmetric coupling between the dots and the respective lead, the current 
noise is sub-Poissonian for resonant tunneling, but super-Poissonian in the vicinity of the resonances. 
Our results indicate that the interaction between different channels together with phonon emission 
and absorption are responsible for the shot noise characteristics. The observed asymmetry of the 
peaks at low temperatures stems from spontaneous emission. 

PACS numbers: 73.63.Kv, 73.23.Hk, 72.70.+m, 73.40.Gk 
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I. INTRODUCTION 

The aim of controlling and manipulating nanoscale de- 
vices requires good knowledge of the processes involved 
in the electronic transport through open quantum sys- 
tems. The increasing success in accessing single elec- 
tron states in semiconductor quantum dots and the un- 
avoidable presence of lattice vibrations in such devices 
obliges one to consider dissipation caused by electron- 
phonon interaction^. The study of the electronic cur- 
rent fluctuations provides further information about the 
system^. E.g., from the investigation of shot noise — a 
consequence of the charge discreteness — we know about 
deviations from Poissonian statistics indicating correla- 
tions between tunneling events. 

A particular example for Poissonian statistics is the 
electron transport through a point contact, for which 
all tunneling events are statistically independent. For 
resonant tunneling through single quantum dots, this is 
no longer the case: As long as an electron populates 
the quantum dot, no further electron can enter and, 
consequently, tunneling events are anti-bunched. How- 
ever, when several of such transport channels conduct in 
parallel and are coupled capacitively, the current noise 
becomes super-Poissonian, as has been demonstrated 
experimentally 5 '-. This means that electrons tend to be 
transferred in bunches, which at first sight is counter- 
intuitive if one thinks in terms of the Pauli exclusion 
principle. The phenomenon can be understood in terms 
of Coulomb interactions between electrons in different 
channels, so that an electron in one channel suppresses 
the transport through the other— and one observes dy- 
namical channel blockade (DCB). Consequently the elec- 
tron transport through one dot occurs in bunches during 
lapses of time when the other dots are empty. 

In a recent experiment— with transport channels that 
consist of double quantum dots (see Fig. [1]), intriguing 
noise properties have been observed: By slightly modify- 
ing the source-drain voltage, the levels of a double quan- 
tum dot can be tuned across a resonance which yields a 
current peak at whose center, the noise is sub-Poissonian. 



In the vicinity of such resonances, by contrast, the noise 
is super-Poissonian such that the Fano factor assumes 
values up to 1.5. This structure becomes washed out 
with increasing temperature, indicating the suspension of 
DCB by the interaction with substrate phonons. In this 
work, we show that a model with a single transport chan- 
nel qualitatively reproduces this behavior. For a quanti- 
tative agreement with the experimentally observed Fano 
factor and temperature dependence, however, the capac- 
itive coupling to a second, almost identical channel is 
found to be essential. 



II. MODEL 

We start out by modelling a single transport channel 
of the setup sketched in Fig. [T] The double quantum 
dot coupled to fermionic leads and substrate phonons is 
described by the Hamiltonia n 1 1 1 12 



H = Hn 



Hi 



cads 



Ht + #e-ph + Hph 



(i) 



where H = J2i=L,R £ i n i + Un h n R - 0(c[cr + c K c L )/2 
describes the coherent dynamics inside the double dot 
and ni denotes the population of dot I = L,R. Hence- 
forth we will assume that the Coulomb repulsion U is so 
strong that only the zero-electron state |0) and the states 




FIG. 1: (color online) Sketch of the transport through two 
parallel double quantum dots measured in Ref™. Both trans- 
port channels are capacitively coupled. The source-drain volt- 
age shifts the relative position of the levels by eV, albeit it is 
so large that all levels lie within the voltage window. 
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with one electron in the left or the right dot, |L) and |R), 
play a role. The leads and the phonons are described 
by fficads = J2i,k £ iknik and H ph = Y^ v h^ v a\a v , respec- 
tively, where nik is the electron number in state k in lead 
I and a„ is the annihilation operator of the isth phonon 
mode. The interaction with the double dot is given by the 
tunneling Hamiltonian Ht = Y2i fc(7'^L Ci + h.c.) an d the 
electron-phonon coupling i/ c -ph = Tl v { n L — riR)X v (al + 
a v )^. By tracing out the leads and the bath within a 
Born-Markov approximation, we obtain for the reduced 
density matrix the equation of motion 
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p = Cp = (C + C T + Cc-ph) p. 



(2) 



Introducing for the density matrix the vector notation 
P = (Poo, Pll, Plr, Prl, Prr) T , the Liouvillian reads 
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where the detuning 5 



J 

(3) 

£r — £l — eV depends on the 
source-drain voltage (or on the gate voltages in lateral 
quantum dots) and E 2 = S 2 + ft 2 . For the phonons, we 
assume an Ohmic spectral density J(u>) = tt ^„ X 2 S(u> — 
lo,j) = 2irau: 1 , so that their influence is determined by 
the coefficients 



A± = 2Tran±2TTaSn( 
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where 7 = Tr/2 stems from the additional decoherence 
associated with the tunneling to the leads. In consistency 
with the experiment of Ref,~, we have assumed that the 
voltage is so large that the Fermi level of the left (right) 
lead is well above (below) the energy of the left (right) 
dot level. Therefore it is sufficient to consider only unidi- 
rectional transport from the left lead to the right lead 13 
described by the effective tunneling rates T; which are 
proportional to \ji\ 2 - 

Within the same approximation, one can derive for the 
current, defined as the time derivative of the charge in 
the right lead, the expression / = etr sys [(j7+ — j7L)po]) 
where po denotes the stationary solution of the master 
equation ^ and J7± are the superoperators describing 
the tunneling of an electron from the right dot to the right 
lead and back, respectively. For unidirectional transport, 
they read J- = and J+p = (T R /h)p RR \0) (0|. 

The noise will be characterized by the variance of the 
transported net charge which at long times grows linear 
in time, (AQp) = St. For its computation, we intro- 
duce the operator tri cads+p h(A r R/9totai)^ i ^ which resem- 
bles the reduced density operator and obeys 



< 



?3 



20 



-20- 




r=o 

T=1.46K 
T=2.92 K 
T=5.83 K 




V [mV] 



FIG. 2: (a) Differential conductance, current (inset), and (b) 
Fano factor through a single channel for various temperatures 
and T L = 0.025, T R = 0.0125, Q = 0.025, e = 0.5 and a = 
0.005 (in meV). The inset of panel (b) shows the Fano factor 
for dissipation strength a = (solid), 10 -3 (dotted), and 10 -2 
(dashed) at zero temperature. 



One can show^ that £ has a divergent component which 
is proportional to po and does not contribute to the zero- 
frequency noise S. Thus S is fully determined by the 
traceless part £j_ = Co — Potr£. In terms of po and 
the zero-frequency noise readsi^ 



S 



e 2 ti sys [2(J+ 



J-K± + (J+ + J-)po]. (7) 



c(t) = cat) + (j+-j-)p(t). 



(6) 



A proper dimcnsionlcss measure for the noise is the Fano 
factor F = S/el which equals unity for a Poisson process, 
while a larger value reflects electron bunching. 



III. TRANSPORT THROUGH A SINGLE 
CHANNEL 

Figure [5^ shows the differential conductance and the 
current for various temperatures as a function of the in- 
ternal bias. In contrast to the dissipationless case (a — 
0) 13 ' 16 , the shape of the curve is no longer Lorentzian but 
exhibits an asymmetry. At higher temperatures, the peak 
becomes broader and more symmetric. This behavior is 
also reflected by the noise. In the absence of dissipa- 
tion, the Fano factor deviates from the Poissonian value 
F = 1: For Tl > Tr (as in the experiment) and a — 0, we 
observe an antiresonant behavior with a dip (F rj 0.5), 
which is accompanied by two maxima with values slightly 
above 1. This double peak structure does not appear if 
Tl < T-r. With increasing dissipation strength a and in- 
creasing temperature, the maxima vanish and the Fano 
factor eventually tends to the Poissonian value F = 1. 
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Although this behavior resembles the experimental 
findings reported in Ref^, there are significant quan- 
titative differences. For the maximal peak value of the 
Fano factor, which is assumed in the dissipationless limit 
a — > for 6 = £7 / \/2, we find the analytic expression 

rp i rx\ -i , (r L - r R ) 2 

p[a ~ ' + 202(r L r R + 2Tl r|) + 8 r*r| • (8) 

It implies F p < 5/4, with the maximum assumed for 
Tr <C This means that for a single channel, 

the theoretical prediction for the maximal Fano factor 
is clearly smaller than the value observed in the exper- 
iment even at finite temperature and in the presence of 
dissipation 10 ; cf. inset in Fig. EJd. Therefore, we must 
conclude that the one-channel model does not fully cap- 
ture the experimentally observed shot noise enhance- 
ment. 
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IV. TRANSPORT THROUGH TWO COUPLED 
CHANNELS 

The natural assumption is now that the shot noise 
must be influenced also by the interaction with a sec- 
ond transport channel; cf. Fig. Q] and ReLi£. Thus, 
we now consider two capacitively coupled channels, so 
that the system Hamiltonian reads Hq = J2ii( £ n n n + 
lEi' v Uu'ii'tiuni'i'), where i = 1,2 labels the different 
transport channels. Note that without interchannel in- 
teraction (Uu'w = for i i'), both channels are statis- 
tically independent. Thus, the behavior observed in the 
one-channel case (see Fig. [2]) is repeated at a different 
voltage, but still the Fano factor cannot exceed the value 
5/4. 

In order to simplify the model, we assume that the in- 
teraction Uu'w is huge whenever i = i' or I = I'. Then, 
the system will accept up to two extra electrons pro- 
vided that they are placed in different stacks and differ- 
ent layer a 8 i 17 . This means that we have to consider the 
following seven states (the ith letter refers to channel i): 
the empty state |00), the one-electron states | L0) , |R0), 
|0L), | OR), and the two-electron states |RL), |LR). We as- 
sume that both dots on the right-hand side couple to the 
same lead, while each channel couples to an individual 
phonon bath. Then we derive for the coupled channels a 
master equation of the form © with a Liouvillian given 
by an 11 x 11 matrix. A closer inspection of this Liou- 
villian reveals that — formally — it can be obtained also in 
the following way: One writes the reduced density op- 
erator of the double channel as a direct product of each 
channel, p = p^ ® p^ 2 \ and the Liouvillian accordingly 
as C = £W + £&, where is the Liouvillian © with 
the parameters replaced by those of channel i. In this 
case, ji = (Tyx +F;r)/2, j ^ i. Finally, one removes 
all lines and columns that contain one of the "forbidden" 
states |LL), |RR). 

For self-assembled quantum dots, a realistic assump- 
tion is that all barriers are almost identical, so that Tl/r 



FIG. 3: (a) Differential conductance and (b) Fano factor for 
two coupled channels for various temperatures and the pa- 
rameters T L = 0.025, T R = 0.0125, n = 0.025, a = 0.005, 
ei = 0.5 and si = 0.75 (in meV). The inset shows the tem- 
perature broadening of the current peak. 



and do not depend on the channel index i. By con- 
trast, for the internal bias £j = £jR — £jL) wc will nn d that 
already small differences play a role, so that we have to 
maintain the channel index i in the effective detunings 
Si = Si — eV. For unidirectional transport, the current 
operators now read J_ = 0, while J + = + acts 
on the reduced density operator as J + p = (Fr/S)[(pro + 
Por)|00)(00| +p RL |0L)(0L| + p LR |L0)(L0|]. 

In the absence of the phonons, we find the scenario 
discussed already in Ref4: The Fano factor exhibits two 
peaks, but their origin is now different than in the one- 
channel case. If both double quantum dots become res- 
onant at different source-drain voltages, an electron in 
the double dot that is out of resonance has only a small 
probability to tunnel through the central barrier. There- 
fore, the non-resonant double dot will mostly be occupied 
with one electron and thereby block the other channel, 
so that the current peaks becomes smaller than in the 
one-channel case; cf. insets of Figs. [5^, and [3^. Whenever 
the nonresonant channel is empty, however, the resonant 
channel will transmit a bunch of electrons, so that even- 
tually the noise is super-Poissonian. 

Figures [3] and [4] show the corresponding current and 
the Fano factor in the presence of dissipation for two dif- 
ferent configurations. We observe two striking features 
which are in accordance with the experimental results 
of Refji -: First, dynamical channel blocking is less pro- 
nounced at higher temperatures and, second, the struc- 
ture of the Fano factor exhibits a clear asymmetry. 

This behavior can be explained within the follow- 
ing picture: Let us consider, for instance, the situation 
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FIG. 4: (a) Differential conductance and (b) Fano factor for 
two coupled channels for various temperatures and the param- 
eters F L = 0.11, T R = 0.055, fi = O.ll, a = 0.005, ei = 0.5 
and £2 = 0.75 (in meV). The inset shows the temperature 
broadening of the current peak. 




FIG. 5: (color online) Phonon-assisted channel opening: 
The blocking electron in the off-resonant channel can tunnel 
through the interdot barrier after phonon absorption (left) or 
emission (right) and thereby open the resonant channel. 



sketched in Fig. [5]where z\ < £2- When the source-drain 
voltage puts the first double quantum dot in resonance, 
i.e., Si = 0, the second double dot is still above resonance 
(62 > 0), thus, blocking the resonant channel. If now the 
electron in channel 2 absorbs a phonon, the blockade is 
lifted. On the other hand, when 82 = 0, double dot 1 one 
is already below resonance (S\ < 0) and phonon emis- 
sion can resolve the blockade. Both processes are more 
frequent the higher the temperature, so that dynamical 
channel blockade is eventually resolved. The fact that 
emission is more likely than absorption, explains the ob- 
served asymmetry and its reduction with increasing tem- 
perature. We emphasize that this does not rely on dif- 
ferences in the interdot hoppings f2i, in contrast to the 
mechanism of RefA 

This phonon-induced channel opening is also mani- 
fested in the enhancement of the current shown in the 
inset of Fig. [3^,. The current peak becomes larger with in- 
creasing temperature and experiences a slight shift away 
in its location. At low temperatures it tends to be around 



the larger resonance voltage {eiI&) which is driven by 
phonon relaxation. As phonon emission becomes impor- 
tant with temperature, the current peak becomes larger 
and shifts towards eVo = (e± +62)/% coinciding with the 
maximal current voltage in the absence of dissipation. 

This effect is weaker when considering stronger tunnel- 
ing couplings: the Fano factor is reduced by DCB lifting 
and the current peak remains centered at eVo when in- 
creasing temperature, which only affects to its boadening, 
cf. Fig. HI 

The observed behavior reproduces rather well the mea- 
surements reported in RefJ^, but there is still one differ- 
ence: In the experiment, the Fano factor far from reso- 
nance is clearly smaller than 1, while for the two-channel 
model, it tends to be Poissonian. This can be explained 
by leakage currents Ik that inevitably flow through the 
whole sample, but have been ignored so far. We assume 
that the leakage currents are statistically independent of 
each other and of the coupled double dots considered. 
Then we can write both the current and the noise of the 
complete sample as a sum of the independent channels: 

Ssys + ^2kFkIk, 



-^sample — -*sys + J2k Ik and S, 



sample 



where is the Fano factor associated to Ik- If the leak- 
age currents stem from resonant tunneling through single 
quantum dots or double dots far from resonance, Fi < 1 
(RefJ^) and, thus, the total Fano factor is decreased: 
-Sample = Sample /Sample < F sys . However, since there 
are about 10 6 leakage channels^, it is not possible to 
estimate their effect more precisely. 



V. CONCLUSIONS 

To summarize, we have studied the effect of electron- 
phonon interaction in the transport through double 
quantum dots systems, predicting super-Poissonian shot 
noise whenever the source-drain voltage tunes a double 
dot close to resonance. The corresponding Fano factor 
exhibits an asymmetric double-peak structure which be- 
comes less pronounced with increasing temperature. In 
order to obtain the experimentally observed values^ for 
the Fano factor, we have assumed that two transport 
channels arc so close that they can block each other. The 
temperature dependence of the double peaks have been 
explained by the suspension of dynamical channel block- 
ing via phonon emission or absorption. We attribute the 
sub-Poissonian noise observed far from resonance to the 
appearance of independent leakage currents. So exper- 
iments with devices where these systems were isolated 
from these additional noise sources are highly desirable. 
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